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Abstract 

The diffusion-controlled reaction kA ^ is known to be strongly dependent on fluc- 
tuations in dimensions d < dc = 2/ {k — 1). We develop a field theoretic renormalization 
group approach to this system which allows explicit calculation of the observables as ex- 
pansions in e^/e^"!), where e = dc — d. For the density it is found that, asymptotically, 
n ~ Akt~'^/'^. The decay exponent is exact to all orders in e, and the amplitude Ak is 
universal, and is calculated to second order in e^/C^"!) for /c = 2,3. The correlation func- 
tion is calculated to first order, along with a long wavelength expansion for the second 
order term. For d = dc we find n ~ Ak(lnt/ty/^'^~^^ with an exact expression for Af^. 
The formalism can be immediately generalized to the reaction kA — * iA, i < k, with the 
consequence that the density exponent is the same, but the amplitude is modified. 
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1. Introduction 



Diffusion controlled chemical reactions are adequately described by mean-field type 
rate equations in higher dimensions, but in lower dimensions the fluctuations become 
relevant . For the reaction kA the critical dimension for fluctuations is conjectured 
to be dc = 2/{k — 1) If d > dc the density n{t) obeys the rate equation 

|n(t) = -^n(^)^ (1) 

with reaction rate constant F. This implies the density will decay asymptotically like 
n ~ {Tt)~^/^'^~^\ For d < dc it is conjectured on the basis of scaling arguments l]^,^, 
rigorous bounds P], and exact results for d = 1 [|6|-|T0i, that n ~ 

t-rf/2_ For d = dc the 

mean- field power law with logarithmic corrections is expected. 

In this paper we apply renormalization group (RG) methods to this system, with 
the goals of verifying the above conjectures and demonstrating universal quantities. The 
formalism developed can be used to calculate the density and correlation function pertur- 
batively in where t = dc — d. This formalism includes infinite sums for each order 

of e^/^^~^\ since the initial density is a relevant parameter and must be summed to all 
orders. 

Previous work in applying RG to this system was carried out by Peliti for the case 
k = 2 [|ll|. Using a field theory formulation of this system, Peliti was able to confirm 
the conjectured decay exponent, and also demonstrate that the reactions A + A ^ ^ and 
A + A ^ A are in the same universality class with regard to the decay exponent and the 
upper critical dimension. Peliti also made the observations that the coupling constants can 
be exactly renormalized to all orders and that there is no wavefunction renormalization 
in the theory. The latter has the consequence that simple scaling arguments can be used 
to extract the decay exponent and the upper critical dimension. However, these scaling 
arguments are not capable of giving other universal quantities in the system, such as 
amplitudes or the asymptotic form of the correlation function. For these one must do the 
complete RG calculation. 

Our formalism enables perturbative calculation of these quantities for general k. For 
example, we find that the density for d < dc is given by n ~ Ak{'Dt)~'^/'^ with 
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and for d = dc 

V(fc-i) + \ i/(fc-i) 



, , , {k-2)\ y'^^-^> flntV 



where V is the usual diffusion constant. 



Recent work in applying RG to this system includes that of Ohtsuki |jT^ , in which the 
density is calculated, although with qualitatively different results than those above. First, 
Ohtsuki predicts that the amplitude for the asymptotic form of the density has the same 
reaction rate constant dependence as the mean-field solution: n ~ for k = 2. Second, 



the leading order term in the e expansion for the density amplitude in [12] is of order unity. 
An RG scheme involving an external source of particles has been developed by Droz and 
Sasvari which leads to scaling functions which confirm the decay exponent. Friedman 
et al. attempted to calculate the density perturbatively, and concluded that it is necessary 
to perform a non-perturbative sum of all orders of no, the initial density, when calculating 
observables ||14|| . This infinite sum is exactly what we do in our calculation scheme. To 
our knowledge there has been no previous satisfactory, complete RG calculation. 

A slightly different field theory formalism for this system was developed in analogy 
with bose condensate calculations |T^,^. This approach leads to a confirmation of the 
decay exponents as well. However, this method is not as readily generalized to an RG 
calculation as is the field theory approach of Peliti. 

The contents of this paper are as follows. In section 2 the system is defined via a 
master equation. This is then mapped to a second quantized representation, and in turn 
to a field theory. In section 3 the renormalization of the field theory and the calculation of 
observables is addressed. The latter requires summing infinite sets of diagrams, for which 
techniques are developed. With the formalism established, the density is then calculated 
in section 4, including correction terms and a discussion of the crossover time scales. An 
alternate method for calculating the leading order amplitude, which does not involve RG, 
is discussed, and its apparent failure in the case k = 2. In section 5 the correlation function 
is calculated, and with it universal numbers for the fiuctuations in particle number, both 
for the total system and for a small volume v. The local fiuctuations in particle number are 
found to be divergent. Also the second moment of the correlation function is calculated, 
giving a correlation length scale. The case d = dc is addressed in section 6, and finally in 
section 7 a summary these results is given, and the generalization to kA iA is discussed. 
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2. The Model 



Consider a model of particles moving diffusively on a hypercubic lattice of size a, and 
having some probability of annihilating whenever k or more particles meet on a lattice site. 
This model is defined by a master equation for P{{n}, t), the probability of particle config- 
uration {n} occurring at time t. Here {n} = (tii, n2, . . . , titv), where is the occupation 
number of the ith lattice site. The appropriate master equation is 



d_ 
di 



P({n}, ^) = ^ 5^ (^e + l)P{. . . , - 1, Tie + 1 . . . , t) - n,P({n}, t) \ 

+Xj2\i^i + k){ni + k-l)... {m + 1)P{. . . , 71, + /c, . . . , t) (5) 
-n^{n^-l)...{ni-k + l)P({n}, t) 



where i is summed over lattice sites, and e is summed over nearest neighbors of i. The 
first curly brackets piece describes diffusion with diffusion constant V, and the second 
annihilation with rate constant A. The P{{n},0) are given by a Poisson distribution for 
random initial conditions with average occupation number uq. 

This master equation can be mapped to a second quantized operator description, 



following a general procedure developed by Doi |17]. To summarize briefiy, operators a 
and are introduced at each lattice site, with commutation relations [ai,aj] = 5ij. The 
vacuum ket is given by ai\0) = 0. The state ket of the system at time t is defined to be 

i</'W) = E^(M,i)n(4no). (6) 

{n} i 

Then the master equation (|]) can be written as 

~^^\m = H\m (7) 

with the non-Hermitian time evolution operator 

^ = - J - a.) - A5^(i - {al)')al (8) 

i e i 

This has the formal solution \(f>{t)) = exp{—Ht)\(l)(0)) . 
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To compute averages it is necessary to introduce the projection state 

N 



(9) 



Then for some observable A{{n}), 



((A(t))) = 5^A(M)P(M,t) 

{n} 



\Aexp{-tH)\(t>{0)) 



(10) 



where A is the second quantized operator analog of A. Note that ( |aj = ( |. Therefore any 
operator A represented in normal ordered form — where all the al have been commuted 
to the left — can be written entirely in terms of the a^. The operator corresponding to 
the density is simply a^, while the correlation function C{xi,Xj) is given by a\aia^-aj 
or aiSij + aiaj. The importance of the S function term will be shown later when the 
renormalized correlation function is calculated. 

The second quantized equation can in turn be mapped to a path integral, with vari- 
ables i/'i, ipi at each lattice site, via the coherent state representation []T^,^. The action 
corresponding to (|^) and (H) is 



- no'i/'i(0) - i/jiit) 



- A(l - 



The last two terms reflect the Poisson initial conditions and the projection state. The path 



integral form of is then 

{{A{t)))=^ I \{di;di^,A{m)e-^^^^^''^- 



(12) 



The normalization constant is given by N'~^ = / Hi c^V'ifiV'ie"'^''^'^'*^- 

Next we take the continuum limit via J d'^x/a'^, ipi — + ip{x)/a'^, 'ipi — * V'(a^), 

^0 noa'^, and '^^{i/Je ~ 4') ^ a'^'V^i/j. The initial density is now uq. The diffusion 
constant exhibits no singular behavior in the renormalization of the theory, so it is absorbed 
into a rescaling of time, giving the action 

-t 



dt<'V(at-v^)^-Ao(i-v')V'' 



- noV'(O) - V(i) 



(13) 
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where Aq = W'^a^''-^'^'^. 

Treating (|TBp as a classical action gives the equations of motion 

dS 



and 



{df - + kXai^^-^i;" - nQd{t') = (14) 



-{df + V')V' + kXoilp'' - 1)^'=-' - S{t' - t) = 0. (15) 



Assuming that and ip are spatially uniform gives the solution ip(t' < t) = 1 and equation 
( P^ becomes 

^7/; = -A;Ao^"+no5(t), (16) 

the mean- field rate equation. It is consistent that the rate constant is /cAq, since Aq 
represents the rate at which the reaction occurs, and the resulting change in particle 
density is proportional to k. Consider shifting -0 by its classical solution: -0 ^ 1 + '0. The 
action which results is (up to an overall constant): 



no^(O) 



(17) 



where Xi = (^)Ao. Note that the boundary terms introduced cancel the '^(t) term in (|T^). 

Averages with respect to this action correspond to physical observables, and are de- 
noted by double brackets. Single brackets are used for averages over the curly bracket part 
of ([T7| ). That is, for some observable A, 

((A(x,t))) = ^A(a:,t)e"°/'''"'^("'°)). (18) 

This is already normalized, since (exp{no'i/'(p = 0)}) = 1. 

The dimensions of the various quantities in (|1^, expressed in terms of momentum, 

are 

[t] = p-^ [^(x)] = / mx)] = / [A,] = (19) 

The couplings become dimensionless at the traditionally accepted value of the critical 
dimension, dc = 2 / {k — 1) [^^. The relative dimensions of i/^ and are arbitrary, but this 
choice is the most natural. Any other choice of dimensions would introduce uq dependence 
into the projection state, and cause the couplings A^ to have different dimensions. 
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3. Renormalization of Observables 

The scheme developed for renormahzing the theory foUows conventional RG analysis 
20|] . In this vein a renormalized coupling is introduced, and shown to have a stable fixed 



point of order e. This is the small parameter of the theory, and not no, which implies 
that the computation of observables requires summing over an infinite set of diagrams, 
corresponding to all powers of no in (p!8|). This infinite sum must be grouped into sets of 
diagrams whose sums give a particular order of the coupling constant. It will be shown 
below that this grouping is given by the number of loops. That is, the infinite set of tree 
diagrams sum to give the leading order term in the coupling, the one-loop diagrams the 
next order term, and so on. However, before addressing the calculation of observables we 
turn to the renormalization of the theory. 

3.1. Renormalization 

To renormalize the theory all that is required is coupling constant renormalization. 
This is because the set of vertices in ( p!7D allow no diagrams which dress the propagator, 
implying there is no wavefunction renormalization. As a consequence the bare propagator 
is the full propagator for the theory. 

To determine which couplings get renormalized one first needs to identify the primi- 
tively divergent vertex functions. A general correlation function with i ip^s and m ip^s has 
the dimension 

[ . . .i^{i)ij{e +l)...i;{e + m))]= (20) 

where (1) = (xi, ti). The Green's function G^^''^\pi, si, . . . s^+m) is calculated by 

Fourier and Laplace transforming the correlation function above, and factoring out overall 
p and s conserving S functions. The dimensions of this quantity are 

^^(£,m)| ^ ^d+2-2^-(d+2)m_ ^21) 

The dimensions of the vertex functions T^^'"^^ are given by the Green's functions with the 
£ + m external propagators stripped off. 

The vertex functions with m < k are those which are primitively divergent for d < dc. 
Since vertices can only connect k tfj^s to some number less than or equal to k i(;^s, then it 
follows that the primitively divergent diagrams have m = k and i < k. 
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A general '0*?/^'^ vertex is renormalized by the set of diagrams shown in fig. |l|. In these 
diagrams the propagator Go{p, t) = {ii'{p, t)i[){—p, 0)) = * for t > 0, Go = for t < 0, 
and is represented by a plain line. Note that this sum is the same for all i, that is all 
vertices renormalize identically. This is a refiection of the fact that there is a only one 
coupling in the theory. These diagrams can be summed to all orders, as noted in [|ll|. In 
(p, t) space the temporally extended vertex function A(p, t2 — ti) is given by 



A(p, t2 - tl) =Xo5{t2 - tl) - Xll{p, t2 - tl 



(23) 



+ Xl / dt'l{p, t2 - t')I{p, t' -tl)-... 

Jtl 

where I{p, t) is the k — 1 loop integral 

/(p,t)=^!/n(|^)(2-)^5(p-x:^'0exp(-x:^'?^)- (24) 

The S function can be written in integral form, which turns the integral into a product of 
k Gaussian integrals. This gives 

I{p, t) = Bkt-^^-^^'^'^e-P^''^ (25) 

where 

Bk^-rj.(-r-] . (26) 



Taking the Laplace transform, A(p, s) = dte~^^X{p, t), makes ( |23|) a geometric sum: 

^ l + XoB,Tie/d,){s + pyk)-^/'^^' ^^^^ 

where the d and k have been exchanged for e and dc- For a general i/j'^ij}^ vertex the Aq in 
the numerator is replaced by A^ = (^)Ao, and the denominator is unchanged. Therefore 
the small s and p form of the vertex function is independent of Aq for all i. 



The vertex function (|27| ) is used to define a renormalized coupling. Using the mo- 
mentum K as a normalization point, we define the dimensionless renormalized coupling to 
be gn = K^^^'^''X{s,p)\s=K.'2,p=o, and the dimensionless bare coupling go = k^^^'^^Xq. The /3 
function is defined by 

(3{gR) = K^gn = -^gn + ^Bkr(-^)gj,. (28) 

OK dr dr. \dr 



It is exactly quadratic in gji and has a fixed point = at 

g*j, = {BkT{e/d,)y\ 



(29) 



The fixed point is of order e. From the definition of (7^, (|27| ), and (^9]) it foUows that 



9r =9o^+9*R \ or 



^0 



9r 



9r 

9R + ^ + --- 



(30) 



1 - 9r/9*r ' 9r 

This wiU be used to exhange an expansion in qq calculated in perturbation theory for an 
expansion in qr. 

3.2. Calculation Scheme 

Next we develop a Callan-Symanzik equation for the theory. Given a correlation 
function 

(31) 



F^^\t, Ao) = (|^(x,t) (^j d''yi;iy,t = 0) 



The lack of dependence on the normalization scale can be expressed via 

F^^\t,K,gR)=0. 



9 , d 



(32) 



From dimensional analysis [F*^"'''] = p'^ implying 



K— 2t— a + ma 

OK ot 



We are interested in the density n{t,nQ, qr, k) = '^^u'q'F^'^^ /m\. Substituting 
(p2D, and summing to get the density gives the equation 



2t^ - dno^ + PigR)^ + d 
ot dno dgR 



n{t,no,gR,K) = 0. 



(33) 



a) into 



(34) 



This is solved by the method of characteristics, and has the solution 



n{t,nQ,gR,K) = {n^t) '^/'^ni^n '^,no{K '^),gR{K ^),k 

with the characteristic equations for the running coupling and initial density 

dfiQ 



(35) 



2t- 



dt 



-dno no{t) = no, 



(36) 



= P{~ga) -g^it) = g^. (37) 
Because of the simple form of the /3 function, the running couphng can be found exactly: 

hoit') = {t/t'f/^no, (38) 

One then sets t' — and plugs the result into (|3^). Notice that in the large t limit 
9r 9r- 

In conventional RG analysis the mechanics developed above is used in the following 
way: one calculates an expansion in powers of g^, and then converts this to an expansion in 
powers of g^i via (^^. As long as the expansion coefficients are non-singular in e, then the 
gR expansion can be related to an e expansion via (]35|). That is, we substitute t —>■ 



no ^ fbQ, gR ^ gji, in the gn expansion, and multiply by the overall factor shown in ( |35|) . 
Then for large t, gn —>■ g'^ giving n{t, tt-q, Aq) as an expansion in powers of e. For a given 
coefficient in the gn expansion we keep only the leading term for large no, since no ~ 
and so the subleading terms in no will correspond to sub- leading terms in t. 

The identification of the leading terms in go is less straightforward than it is in conven- 
tional RG calculations, since the sum over all powers of no must be taken into account. For 
the density, tree diagrams are of order g^n]^^^^ ^'^ for integer i. Diagrams with j loops are 
of order g^n^^^^ . Since the addition of loops makes the power of go higher relative 
to the power of no, we hypothesize that the number of loops will serve as an indicator of 
the order of ^fo- This will be shown to be the case via explicit calculation. 

3.3. Tree Diagrams 

To calculate all possible diagrams of a given number of loops it is necessary to develop 
two tree-level quantities: the classical density and the classical response function. The 
term classical means averaged with respect to the classical action, which is the action (p!7|), 
but with only the ip'tp'^ vertex. The classical density is given by sum of all tree diagrams 
which terminate with a single propagator, as shown in fig. |^, and is represented graphically 
by a dashed line. These diagrams are evaluated in momentum space. From (^8]) it follows 
that the ipit = 0) in the initial state all have p = 0, so all diagrams at tree level have p = 0. 

Shown also in fig. ^ is an exact graphical relation for the infinite sum, which is equiv- 
alent to the mean-field rate equation (|1^) . This can be seen by considering the diagram in 
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position space, and acting with {dt — V^), the inverse of the Green's function Gq, on either 
side of the diagramatic equation. Note that the combinatoric factors involved in attaching 
the fuU density hues to vertices is different than for propagators, which is discussed in 
appendix A. This equation has the exact and asymptotic (large t) solutions 

The asymptotic solution depends on the coupling strength, but not the initial density. 
The response function is defined by 

G{p,t2,h) = m-p,t2)^(^{p,h))), (41) 

and the classical response function is the above quantity with only tree diagrams included 
in the averaging. It is represented graphically by a heavy line, and is given by the sum of 
diagrams as shown in fig. |^. Note that the only p-dependence is that of the bare propagator. 
That is, the density lines all carry no momentum. The time dependence of the propagators 
connecting the vertices cancels to leave only overall dependence on ti, ^2- The vertices are 
now symmetric under interchange, so we can trade the requirement that they be ordered 
for a factor of l/n„! where is the number of vertices. The sum of diagrams is then 
identified as the Taylor expansion of an exponential, giving 

G^iipMM) =e-p'(*^-*i)exp|-A;2Ao 1' dtn,i{tf-^ 

k/(k-i) (42) 



1 

k-i 



/ l + fc(fc-l)< 'Xoh 

The extra factor of k associated with each —kXo vertex is a consequence of the com- 
binatorics (see Appendix A). From ( [T8|) it follows that (('i/'(t)'i/'(0))) = d{{'^{t))) / driQ or 
G{j> = 0,t, 0) = dn{t)/dno. This relation should also hold for the classical density and 
response function, as is the case for the solution above. 

4. Density Calculation 

With the classical or tree-level solutions of the previous section, and the renormaliza- 
tion scheme developed above, the asymptotic form of the density can now be calculated. 
The solution for the tree diagrams in terms of go, or Aq, is given by (P0|). To leading order 
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in gn one just replaces Aq with Qrk^'^^'^'' . For large t the running coupling g^i — > (7^, which 
gives 

nWm = ^ (43) 

The superscript on the density refers to the number of loops in the calculation. The 
asymptotic form of this expression is 



27r(A;- l)A;V(fc-i)f 

The term in parentheses is the leading order term in A^, the amplitude of the t"*^/^ 
component of the density. 

Amplitude Corrections for k — 2 

Next the corrections from the higher loop diagrams are calculated. It will be shown 
that adding a loop makes the sum of diagrams an order g]{^'^ higher. At /c — 1 loops the 
diagrams will contain a singularity in e, caused by the appearance of the first primitively 
divergent diagram. However this singularity is cancelled when the gj^ correction to go in 



( pOD is included in the tree diagram sum. In general the higher order terms in (|30| ) will 
cancel all divergences in the coefficients of the gR expansion. This will be illustrated in 
the one-loop corrections for k = 2. 

The infinite sum of all one-loop diagrams can be written in terms of the classical 
response function found above. The sum of diagrams is shown in fig. ^ Expressing this 
graph in integral form 

n^^\t,nQ, go, k) = 

r (jd (45) 

2 j dt2dtij^Gci{0,t,t2){-2Xo)Gci{p,t2,tif{-Xo)nci{tif, 

where the time integrals are over < ti < t2 < ^- Taking the large no limit of ( ^5]) to 
extract the asymptotic part gives 

nW(t,no,^7o,«) = dt2dt^-^t-^e-^p"^'--'^hl + 0{n-'). (46) 

Notice that this is independent of go, consistent with the prediction that the one- loop 
diagrams are of order (7^ and provide a correction to the leading term in (^1|) . The integral 
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can be done exactly. Expressing the leading piece in terms of g}^, and the rest as an 
expansion in e: 

n(^)(t, no, 90. n) = t"^/^ - ^ + 0(e)) , (47) 

where C is Euler's constant. The correction to the tree-level component due to the sub- 
leading term in go{gR) is 

n^'\t,gn) = ^t"^/^ - ^^t-'^' + 0{gn). (48) 
^9r ^9r 

The singular parts of the g^^ coefficient cancel as advertised. Combining (^T]) and (^SD and 



making use of the Callan-Symanzik solution (|3^) gives 



1 21n87r - 5 , 

The two-loop diagrams are also shown in fig. ^. They all contribute to order g]^. 
Unfortunately we are unable to evaluate diagrams (f,g) due to the complicated time de- 
pendence of the vertices, which prohibits calculation of the 0(e) term in A2. The most 
singular of the diagrams, (b-d), are of order e~^. These diagrams can be calculated and 
the singular pieces cancel as expected. 

Note that the asymptotic, or large no, limits of the classical density and the classical 
response function are of order ng, which implies that the asymptotic time dependence of 
the density, calculated to any number of loops, will be t~'^l'^ . Therefore the decay exponent 
is exact to all orders in e. 

The cancellation of the singularities which appear in the g^ expansion can be most 
easily understood by viewing the correction terms in ( |30| ) as counterterms introduced to 
cancel primitive divergences. That is, considering (Jgr^ = g\l g\ + ^((7^), and calculating 
the first order term in bg^ at tree level gives a diagram similar to fig. § (a), but with the 
counterterm in place of the loop. This diagram, when added to the one-loop diagram, 
cancels the singularity in the g\ coefficient. Two-loop diagrams (b-f) can be viewed as 
primitively divergent loops added to the one- loop diagram (a). The order bg^, terms in 
the one-loop diagram are equivalent to diagrams (b-f) with a counterterm in place of the 
additional loop, and will cancel the divergences in these diagrams. Diagram (g) differs 
in that it is not a primitively divergent loop 'added on' to diagram (a), but it is also 
non-singular. 
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4-2. Amplitude Corrections for k = 3 

The one- and two-loop diagrams for k = 3 are shown in fig. The one-loop diagram 
contains no singularity, and gives the order correction to (0). The asymptotic piece is 
given by the integral 

n^'Ht,no,go,K) = ^/ dt^dt^^t^^-'^'^'^-'^kl/' + 0{n^'). (50) 



Performing the integral and using (|3^) we find the amplitude 



1 /2 

The two- loop diagrams are of order , although, similar to the case of /c = 2, we 
are unable to calculate diagrams (f-i). The only diagram with a singularity is (j) which 

1 /2 

can be calculated to demonstrate that the g^^ coefficient is non-singular as expected. 

4-3. Dressed Tree Calculation 

There exists an alternate method for calculating the leading order amplitude of the 
density which does not require using the RG formalism. However, there is a discrepancy 
between this method, the dressed tree sum, and the RG in the case k = 2. We present 
the dressed tree calculation below, and an explanation for why we believe the RG to be 
correct for k = 2. 

Consider summing the most divergent diagrams for each power of Aq and no. This 
is equivalent to summing the dressed tree diagrams, which are tree diagrams with all 
the vertices replaced by the temporally extended vertex function (^3|). The sum of these 
diagrams, ndt{t), satisfies the diagramatic equation shown in fig. |^, where Udt is represented 
by a dotted line. As with the tree diagram sum, acting on this equation with the propagator 
inverse {dt — V^) gives a differential equation 

dtUdtit) = no5{t) - k I dt'X{p = 0,t-t')ndt{t')\ (52) 

Jo 

Laplace transforming the equation gives 

sn{s) - no = -kX{0,s)n'^{s), (53) 
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where n(s) = dte~^^n{t) and n^{s) = dte~^^n{t)^ . The transform of the vertex 
function A(0, s) is known exactly, and is given by However, the equation is not 

algebraic in n(s), making it difficult to obtain an exact solution. To proceed, we assume 
Udt ~ it"", so that for smaU s, n{s) ~ AV{1 - Also, n^{s) ~ i^r(l - A;a)s^"-\ 

and A(0, s) ~ s^^'^" /{BkT{e/dc)). The transform of Udtit)^ is calculated by imposing a small 
t regulator, which is justified as the transform of the exact solution does exist, and then 
taking the small s limit. The amplitude which results is independent of the regulator. 
Substituting these in to (|5^ ) and taking the small s limit of the eqation gives a — d/2, 
and the amplitude 



,k-i - \k-i ^2Mfc-i "2" 



For k ^ 2 the non-singular F functions cancel to leading order in e, with the result A = 
Ak+0{e*^). However, for /c = 2 all the F functions are singular, which has the consequence 
that A2 = 2A2 + 0(e°). In light of this, it seems necessary to find an explanation why this 
particular set of diagrams sums to give the proper leading order term for general k, but 
not for k = 2, if indeed the RG is giving the correct leading order term. 

Consider the set of dressed one-loop diagrams. That is, the set of diagrams given in 
fig. 1^ (a) and fig. |] (a), but again with each vertex replaced by the temporally extended 
vertex function. While it would be difficult to calculate this sum, it is possible to see a 
property specific to /c = 2 that they have. The analog of the classical densities in these 
diagrams is the dressed tree density Udt oc t"*^/^. Therefore for general k there is a time 
integral over t~^'^^'^, or f-'^/i'^-^)-'^'^/'^^ This time integral will be in the form of a Laplace 
convolution integral, similar to ( |5^ . Using a regulated transform as before, the amplitude 
of the small s limit will be proportional to F((/c — 2)/(k — 1) + ke/2). For k ^ 2 this is 
non-singular at e = 0, but for /c = 2 it is of order e~^. Therefore these diagrams are part of 
the leading order amplitude for /c = 2. As a result, it would appear that the discrepancy 
is a consequence of the failure of the dressed tree method, and not of the RG. 

4.4- Crossovers 

There are two crossover time scales in this system, one given by no and one by Aq. 



For the coupling constant crossover we consider the large t expansion of (|39| ) 



-gn = 9*r{i- Xoh-^/"^ + 0(1-'^/"^)) . (55) 
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Including the correction term in the density calculation will generate a Aq dependent 
term proportional to i-d/'^-e/dc _ From (|55|) it follows that the characteristic crossover 
time is given by ~ (e/Ao)'^''^'^- In terms of the constants in the master equation, 
t\ ~ a?''D~^{e'D / a^XY"/'^ . For small e, or large Aq, the time required to reach the fluctuation 
dominated regime becomes small. 

The no crossover is calculated by keeping the order tiq^ terms in the integrals per- 
formed above. These terms will pick up an extra factor of when put into (|3^), so 
the exponent of the leading no dependent term in the density is t~'^. The characteristic 
crossover time is only weakly e dependent, and is given by ~ 'D~^nQ ^^'^ = a^P~^nQ 

If the no crossover occurs first, then for intermediate times ^ t <^ one would 
expect the system to obey the asymptotic form of the mean-field solution. That is, n ~ 
[k{k — 1) Xot]~^ / ^'^~^K If the Ao crossover occurs first it is less clear what the behavior in the 
intermediate regime will be. The contribution from the tree diagrams will be exactly (^), 
which does not become a power law until the no crossover is reached. This is complicated 
even further by the higher order diagrams. 

5. Correlation Function Calculation 

The density correlation function is given by 

C{x,t) = {{{i;{x,t) + S'{x))i;{0,t))). (56) 

where the S function is a consequence of the second quantized operators developed in 
section 2. A Callan-Symanzik equation for the correlation function can be developed in a 
similar fashion as before. Consider the function 



n(m) 



{p,t,Xo)^ Jd''xe-'P-^(^{i;{x,t)+5''{x))i:{0,t)(^Jd''y^j{y,t = 0)^ ^ (57) 



Dimensional analysis gives [F^"^)] = p'^ The correlation function C{p,t) is given by 
nQ^F("^)/m!. This leads to the equation 



' d 9 , d , d 

ot op OUq OQu 



C{p,t,no,gR,K) = 0, (58) 



which has the solution 



C(p, t, no, gR, k) 



= {Kh)-''/'c{p{K-'),t = K-^M^-')r9R{^-').^), (59) 
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with Qfi and no given by ( |38|) and (^9|), and 

Pit') = P\ll- (60) 

Again the calculation of the right hand side of is divided into the number of loops. 
First the connected and disconnected pieces are separated 

C(p, t) = n{t) + g{p, t) + d\p)n{tf. (61) 



The first term on the right hand side is a consequence of the 5 function in (pq), and is 
considered part of the connected correlation function. The disconnected tree-level graphs 
are of the order QqUq^^^'^ ^\ and represent the leading order terms in the correlation 
function. This is reasonable, as the classical solution of this system corresponds to the 
absence of correlations. The connected tree-level diagrams, which are the leading terms in 
g{p, t), are of order QqUq^^^'^ ^\ and represent the leading corrections due to fluctuations. 
The tree- level and one- loop diagrams for g{p, t) in the case k = 2 are shown in fig. ^ 
Diagram (a) can be calculated explicitly to give the leading term 

9{p,t) = --^t-'^/'f2{p't)+0{e') 

I (62) 

J L ^ 

^^^""^ ~ 4x3 ^ 4x3 2x2 ^ 2x' 

The function /2(x) is regular at x = 0, with /2(0) = 1/3. For large x, /2(x) ~ l/(2x). 

We are unable to evaluate the one-loop diagrams analytically for general p, although 
it is possible to calculate an expansion in p^ , which we have done to order p^ . For the 
connected correlation function, C{j>, t) = n{t) + g{p, t). 



C{p, t) 



1 91n87r-26 / 1 151n87r-19 



6ne IOStt \24ne 720n 



+ — TT^T^ + 1 + ^—^^ — 1 P^t + 



t-'^/2 + 0(e). (63) 



With the expansion above it is possible to calculate the second moment of C'(x, t), giving 
a length scale for the correlations. For C{p,t) = A + Bp^ -f . . . the second moment 
—^2 = j d'^x x^C{x,t) / J d^xC{x,t) = —2B/A. The negative sign in the definition of ^ is 
required since the second moment is negative, indicating that the particles are negatively 
correlated at larger distances. For k — 2 the length ^ is given by 
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The correlation function can be used to calculate the fluctuations in the density. 
For example, the fluctuations in the local density are given by integrating C{p,t) over p. 
However, the p-independent term causes this integral to diverge. One can consider the 
fluctuations of the average particle number of flducial volume v. This is given by 

{SNyf = v f dxC{x,t) =vn{t) +0{v^), (65) 

Jv 

where translational invariance is assumed. The order v contribution originates from the 
S function in (^6[). For small v the fluctuations go as SNy ~ a/ vn{t), which is universal. 



Also, 5Ny/Ny ~ l/^Jvn{t)^ which diverges as v goes to zero, consistent with the local 
fluctuations being divergent. 

The fluctuations in the total number of particles is given by VC{p = 0, t) where V is 
the volume of the system. When divided by the square of the average number of particles, 
y^n(t)^, this gives 

(SNV^^ /Stt 367rln87r - 767r 2 h/2 , n 

-y= — e e2 + C>(e^) (66) 



A^2 \^ 3 6 

Note that all these fluctuation terms would be negative if the b function term were ne- 
glected. That is, {{'i/j{x)'^)) < 0, a demonstration that the flelds introduced via the path 
integral formulation of [|18[ are complex. 

The diagrams contributing to g{p, t) for k = 3 are shown in flg. ^. The leading order 
term for the connected part is 

1/2 



12ne 



e-2^3v^erfl(y2^) 3 1 



(67) 



where erfl(a;) = — zerf(za;) = (2/y^) dye^ . The function f'i{x) is also regular, with 
/3(0) = 2/5 and /3(x) ~ l/{2x) for large x. 

The one-loop diagrams can be calculated as an expansion in p^, with the net result 

^/ 3 / VsV^^ 81v^ 288 [2 

^^^'')= lols^J +^60^ + 875V^ 



.1,2 ^ . (68) 

, V3 \ 9v^ 27546 [2 , ^ 
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In this case the sign of the second moment of the correlation function depends on e. For 
e < 0.4 the second moment is negative, and the resulting length scale is given by 

1/2 



/ 4 2^/?,\ nUn 2217\ ./o , 



(69) 



The fluctuations in total particle number are given by 



iv2 \^ 5 ^ ^ V 875 400 y 7 

6. d = dc 

In general, when d < d^ certain relevant parameters determine the critical exponents 
of the system. When d = dc these parameters become marginally irrelevant. In such a 
case the exponents are given by mean-fleld theory, but with logarithmic corrections. In 
our system the marginally irrelevant parameter is the coupling Aq. 



When d = dc the Callan-Symanzik solution (|35|) still holds, although with a different 
running coupling. The (3 function can be calculated either with a cutoff which is taken to 
infinity or by taking e ^ in ( PBD with the same result: (3{gii) = 2Bkg%- This gives the 
running coupling 

~9Ri^~") = 1^ ^rf,. 2,v (71) 
1 + guBk ln{KH) 

For large t the coupling goes to zero, which is the only fixed point of the (3 function. Using 
the asymptotic form gji ~ {Bk ln(K^t)}~^ in the tree-level sum gives 

"(*'-(4w) (¥) ho((,„r (72) 

Higher order terms in gn will give sub-leading time dependence, so this represents the full 
leading order amplitude. Notice that the correction terms are only an order ([nt)~^/^'^~^^ 
smaller, which will make time required to reach the asymptotic regime large. 

The same procedure gives an exact expression for the leading term in the correlation 
function as well. For k = 2 

Cip^t) = ^[l-hiph)) (i^) [l + 0((lnt)-i)] (73) 



and for k 



CiP^i)=i^) ' {^-fsip't)) (^)'^'[l + 0((lnt)-V2)]. (74) 
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7. Summary and Generalization to kA £A 

With the RG calculation developed above we are able to calculate various universal 
quantities for this system. These include the amplitude of the asymptotic density for 
d < dc, given by (^9|), (|51]), and ([72D, and the connected correlation function. Also universal 
are the fluctuations in total particle number and the fluctuations in particle number in a 
small volume v. 

The density amplitude for = 2 can be compared to the exact solution ior d = 1 
of A2 = (87r)-i/2 ^ 0.20 0]. Putting e = 1 in our expansion yields A2 = 0.08 + 0.03 + 
. . .. The agreement is less than satisfactory, indicating that the e expansion will not be 
quantitatively accurate to e = 1. However, the e expansion provides the only systematic 
derivation of universality and scaling. 

Our results can be immediately generalized to a coagulation reaction kA iA, £ < k. 
The only change in the field theory is the vertices in (|T7|): 



Ao©-Ao(^), i<i 

(75) 

Ao(J), i>t 

The renormalization follows identically. For example, the leading term in the amplitude, 
given by (|4l), is generalized to 

^M=(^j Ak + 0{e'). (76) 

This proportionality is not generally true for all terms in the e expansion, although it does 
happen to hold when k = 2. To see this consider a rescaling —>■ 6-0, — > t(j/b, and 
riQ —>■ briQ in the action ( |T7D . The only terms changed by such a rescaling are the couplings 
Xi V'~'^\i, which for /c = 2 is only the coupling Ai. Starting from the theory A + A ^ 
and making the scale transformation with 6 = 2 gives exactly the theory for A + A ^ A. 
As a consequence, the density for A + A ^ A, starting from an intial density of no, will for 
all times be exactly twice the density of the system A + A ^ with initial density of 7io/2. 
This result agrees with the recent exact solution of a particular model of A + A ^ (0, A) in 
d = 1 [0, although it should be noted that this relation is not truly universal for all times. 



as it only holds when the irrelevant couplings are excluded. The asymptotic amplitude is 
universal, and so the relation ^2,1 = 2^2,0 is exact to all orders in e, and independent of 
the initial densities. 
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For k = 3 such a simple relation does not hold. We can consider all three theories, 
i = 0,1,2, combined with relative strengths ro,ri,r2, where ~ l- '^^^ rescaling 

defined above will relate two systems with different in that the densities will be identical 
up to a rescaling. However this rescaling only removes one degree of freedom from the two 
independent variables, so unlike k = 2, one cannot necessarily scale one theory into another. 
Considering tq and ri, we find 

ro{b) = (1 - bf + b{2b - l)fo + b{l - b)fi (77) 

n (b) = (3 - 6) (6 - 1) + 26(1 - 6)fo + 6Vi (78) 

where fg, fi are the values of ro, ri prior to rescaling. Consider the system which is purely 
£ = 0, or fo = 1, fi = f2 = 0. For any b 1 then ri(6) < 0, which implies that there is 
no combination of systems with different i which is equivalent to £ = up to a rescaling. 
This is not the case for the pure £ = 1 system. This system can be rescaled from 6 = 1 to 
b = 3/4. At the latter point one has tq = 1/4, ri = 0, and r2 = 3/4, so this combination 
of systems, with an initial density of 3no/4, will give exactly 3/4 the density of the £ = 1 
system at all times. Similarly, starting with f2 = 1 the system can be rescaled from 6=1 
to 6 = 3/2. At the latter point ro = 1/4, ri = 3/4, and r2 = 0. 

It should be noted that the correlation function will not be identical up to a rescaling 
for any of the systems described above. This is a consequence of the fact that the correlation 
function contains both t/j and i/j^ pieces. 

While the reaction considered here is not as generally interesting as that of A + B — > 0, 
it is a suitable starting point for developing the application of RG methods to these systems. 
A similar approach may be applicable to the reaction mA + nB — 0, a system where the 
universality classes appear to depend on the nature of the initial conditions [pT]-p3 . 
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Appendix A. Symmetry Factors 

Diagrams which contain the classical density or the classical response function are 
representations of infinite sums of diagrams. While they resemble ordinary perturbation 
theory diagrams, they differ in combinatorics. When calculating the Wick contraction com- 
binatorics one treats propagators as distinguishable, although the resulting combinatoric 
factor is then cancelled by a factor which is absorbed into the definition of the coupling 
constant. Our diagrams differ from this in two ways. First, the classical density is attached 
to vertices as an indistinguishable object. This will be demonstrated below. Second, we 
have chosen to introduce in the coupling constants no pre-adjusted combinatoric factor. 
This is merely a matter of convention, and is motivated by the indistinguishability men- 
tioned above, and by the direct relation of the coupling constant to the parameters used 
in the master equation. 

The indistinguishability of the density lines can be demonstrated by considering the 
contraction of k i/^'s, representing a vertex, with the infinite sum which is the initial state. 

m=0 

oo k . . (A.l) 

m=0 ' mi,...,mk i=l V / 

mi + ...+mk=m 

where C*^^ = m!/(mi! . . .m/j!) is the number of ways to partition m objects into k 
distinct boxes. The sums can be replaced with unrestricted sums over mi . . .ruk, and the 
above expression factors completely, giving 

u'))ci = m'ci (A.2) 

The significance of ( |A.2|) is that there is no k\ prefactor. The k classical density lines which 
are connected to the vertex are effectively indistinguishable. 

In calculating the classical response function it is necessary to consider attaching one 
propagator and k — 1 density lines to a "(p^ vertex. This brings in a factor of /c, for the 
number of distinguishable ways the propagator can be attached. The remaining k — 1 
densities follow through the same combinatorics as that shown above, and contribute a 
factor of 1. 

In general, where the classical response function appears in a diagram it can be treated 
as a propagator for combinatorics. The exception to this situation is in diagrams such as 
fig. ^ and fig. 1^, diagrams (d). Here the symmetry of the two disconnected branches will 
result in the branches attaching as indistinguishable objects. 
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Figure Captions 



Fig. 1. Sum of all the diagrams which contribute to X{p, t2 — ti). Shown here is the case 
k = 3, i = 1. These diagrams can be summed exactly, and are the same for all i. 

Fig. 2. The classical density, represented as a dashed line, is given by (a) the complete 
sum of tree diagrams, and (b) an integral equation. The latter is equivalent to 
the mean- field rate equation. Shown here is the case k = 2. 

Fig. 3. The response function, shown as a heavy line, is given as a sum of the bare 
propagator plus a term with a single vertex connecting k — 1 full density lines, 
plus a term with two vertices, and so on. Shown here is A; = 3. These diagrams 
can be summed exactly. 

Fig. 4. One- and two-loop diagrams for k = 2. By using the response function all such 
diagrams are included. Diagram (a) is used to calculate the amplitude correction. 

Fig. 5. One- and two-loop diagrams for A; = 3. Diagram (a) contains no e singularity, 
and is used to calculate the amplitude correction. 

Fig. 6. Exact diagramatic equation for ndt(t), the sum of the dressed tree diagrams. 

Fig. 7. The diagrams for the connected correlation function at tree level and one loop, 
for k = 2. 

Fig. 8. The diagrams for the connected correlation function at tree level and one loop, 
for k = 3. 
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